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v^=10.52 GeV 

Abstract 

Using the CLEO detector at the Cornell Electron Storage Ring, we have 
made a measurement of R= '^^'^g+g^^^t''^-^^'* = 3.56±0.01±0.07 at -y/s=10.52 
GeV. This implies a value for the strong coupling constant of (10.52 GeV) = 
0.20 ± 0.01 ± 0.06, or a^(Mz) = 0.13 ± 0.005 ± 0.03. 
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I. INTRODUCTION 



The measurement of the hadronic production cross section in e^e^ annihilation is perhaps the 
most fundamental experimentally accessible quantity in quantum chromodynamics (QCD) due to 
its insensitivity to the fragmentation process. The measured hadronic cross-section is generally 
expressed in terms of its ratio R to the point cross-section for fJ-^fJ-^ production. In QCD, R is 
directly proportional to the number of colors, depends on quark charges, and varies with energy, 
both discreetly as quark mass thresholds are crossed, and gradually as the strong coupling constant 
as "runs". So, historically, R measurements have been valuable in verifying quark thresholds, 
charges, color-counting, and the existence of the gluon. 

The theoretical prediction for R is 

R = ii(o)(l + as/TT + C2(a,/7r)2 + Csias/nf). (1) 

A calculation appropriate at LEP energies obtained C2 = 1.411 and C3 = —12.68 |j^] for five active 
flavors, in the limit of massless quarks. A recent calculation, applicable to the T mass region (four 
active flavors), has included corrections due to the effects of quark masses and QED radiation 
to obtain 6*2=1. 5245 and C3=-11.52 at y^=10 GeV p. The effect of including these additional 
corrections is a difference of approximately 0.3% in the prediction for R at this energy. i?(o) is 
the lowest-order prediction for this ratio, given by i2(o) = ^c^iQi, where Nc is the number of 
quark colors; the sum runs over the kinematically allowed quark flavors. Just below the T(4S) 
resonance, where bb production is kinematically forbidden, the lowest-order prediction is therefore 
obtained by summing over udcs quarks, yielding i?(Q)=10/3. An experimental measurement of R 
can therefore be used to deduce a value for as- In this Article we present a measurement of R using 
the CLEG detector operating at the Cornell Electron Storage Ring (CESR) at a center-of-mass 
energy ^/s=W.52 GeV. 

II. APPARATUS AND EVENT SELECTION 

The CLEG II detector is a general purpose solenoidal magnet spectrometer and calorimeter 
IQ]. The detector was designed to trigger efficiently on two-photon, tau-pair, and hadronic events. 
As a result, although hadronic event reconstruction efficiencies are high, lower-multiplicity non- 
hadronic backgrounds require careful consideration in this analysis. Good background rejection is 
afforded by the high precision electromagnetic calorimetry and excellent charged particle tracking 
capabilities. Charged particle momenta are measured with three nested coaxial drift chambers 
with 6, 10, and 51 layers, respectively. These chambers fill the volume from r=3 cm to r=100 cm, 
where r is the radial coordinate relative to the beam (z) axis, and have good efficiency for charged 
particle tracking for polar angles |cos^| <0.94, with 9 measured relative to the positron beam 
direction {+z). This system achieves a momentum resolution of (Sp/p)'^ = (O.OOlSp)^ -|- (0.005)^, 
where p is momentum in GeV/c. Pulse height measurements in the main drift chamber provide 
specific ionization resolution of 6.5% for Bhabha events, giving good K/tt separation for tracks with 
momenta up to 700 MeV/c and approximately 2 standard deviation resolution in the relativistic 
rise region. Gutside the central tracking chambers are plastic scintillation counters that are used as 
fast elements in the trigger system and also provide particle identification information from time-of- 
flight measurements. Beyond the time-of-flight system is the electromagnetic calorimeter, consisting 
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of 7800 thallium-doped cesium iodide crystals. The central "barrel" region of the calorimeter covers 
about 75% of the solid angle and has an energy resolution of about 4% at 100 MeV and 1.2% at 
5 GeV. Two endcap regions of the crystal calorimeter extend solid angle coverage to about 98% of 
47r, although with somewhat worse energy resolution than the barrel region. The tracking system, 
time-of-flight counters, and calorimeter are all contained within a 1.5 T superconducting coil. 

To suppress rr, 77, low- multiplicity QED, and other backgrounds while maintaining rela- 
tively high qq event reconstruction efficiency, the following requirements are imposed to define our 
hadronic event sample: 

1. At least 5 detected, good quality, charged tracks (Nchrg) >5. 

2. The total visible energy Evis {=Echrg + Eneutrai) should be greater than the single beam 
energy: Evis > Ebeam- 

I pmiss I 

3. The z-component of the missing momentum must satisfy: ' ^ ^ ' <0.3. 

4. To suppress events originating as collisions of beam particles with gas or the vacuum 
chamber walls, we require that the reconstructed event vertex (defined as z^rtx) be within 
5.5 cm in z {z defined above as the e"*" beam direction) and 2 cm in cylindrical radius of the 
nominal interaction point. 

5. In addition to these primary requirements, additional criteria are imposed to remove back- 
grounds remaining at the ~1% level, as well as to suppress events with hard initial state 
radiation, for which theoretical uncertainties are large. These are: 

(a) No more than two identified electrons are in the event. 

(b) The ratio R2 of the second to the zeroth Fox- Wolfram moments ||5| for the event should 
satisfy R2 <0.9. 

(c) The ratio of calorimeter energy contained in showers that match to charged particles 
relative to the scalar sum of the momenta of all the charged particles in the drift chamber 
(^) must be less than 0.9. 

(d) The most energetic photon candidate detected in the event must have a measured energy 
less than 0.75 of the beam energy (x^ = — <0.75). This requirement reduces the 
uncertainty from radiative corrections. 

Figures |^-^ show comparisons between our candidate hadronic event sample in data and the Monte 
Carlo sample (arrows indicate our cut values). Simulated hadronic events are produced using the 
JETSET 7.3 qq event generator run through a full GEANT-based detector simulation. Tau- 
pair events use the KORALB event generator in conjunction with the same detector simulation. 
In all of these comparison plots, both the data and the 'Monte Carlo sum' have been normalized to 
unit area in the 'good' acceptance region. We also use this Monte Carlo event sample to determine 
the efficiency for qq events to pass our hadronic event selection requirements. We note that, 
according to the Monte Carlo simulation, the trigger inefficiency with the default event selection 
criteria is less than 0.1%. This has been checked with the data by counting how many hadronic 
events triggered only a minimum bias, pre-scaled trigger line. 
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FIG. 1. Normalized visible energy distribution for data (solid), qq Monte Carlo (dashed), and 
rr Monte Carlo (dotted). Sum of qq plus rr Monte Carlo is shown as crosses. 
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FIG. 2. Normalized charged multiplicity distribution for data vs. Monte Carlo. 
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FIG. 3. Distribution of z-component (i.e., direction cosine) of missing momentum p'^'^^^ j^p^' 
for data vs. Monte Carlo (this variable is not cut on). 
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FIG. 4. Ratio of P^^^^ j EmsiUe for data vs. Monte Carlo. Two-photon collisions, and beam-gas 
interactions tend to populate the regions away from zero and towards ±1 in this plot. 
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FIG. 5. Z-component of thrust axis for data vs. Monte Carlo (this variable is not cut on). 
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FIG. 6. Ratio of transverse momentum relative to visible energy (this variable is not cut on). 
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FIG. 7. Comparison of data vs. Monte Carlo spectrum of most energetic photon observed in 
event. 

III. BACKGROUNDS 

After imposition of the above event selection criteria, we are left with a sample of 4.00 x 10^ 
candidate hadronic events. Small backgrounds still remain. These are enumerated as follows: 

1. Backgrounds from e'^e" t+t^(7) events are subtracted statistically using a large Monte 
Carlo sample of KORALB tau-pair events. These events comprise (1.3 it 0.1)% (statistical 
error only) of the sample passing the above event selection criteria. 

2. Backgrounds from the narrow T resonances can be explicitly determined from data us- 
ing e+e- ^ 7T(35/25); T(35/2S) ^ 7r+7r-T(15); T{1S) l+l- events. These events 
are distinctive by their characteristic topology of two low-momentum pions accompanied 
by two very high momentum, back-to-back leptons; the photon generally escapes unde- 
tected along the beam axis. As shown in Figure |8|, we observe these events as dis- 
tinct peaks in the mass distribution recoiling against two low-momentum pions in events 
also containing two high-energy muons. (The recoil mass is calculated from: MrecoU = 

{2Eheam — -E'tti — Ej^^Y — {p-jr-^ + ^712)^, and therefore neglects the four-momentum of the 
initial state radiation photon.) Backgrounds from QED processes {e~^e~ — > j£~^i~; 7 
e^e~) can be suppressed by requiring that the candidate dipion system not be colinear with 
either of the final state leptons. Knowing the branching fractions Q for T(2S)— > 7r7rT(lS) 
(18.5±0.8%) andT(3S)^ 7r7rT(lS) (4.5±0.2%), the leptonic branching fraction for the T(1S) 
(2.5 lb 0.1)%, and the reconstruction efficiency for such events (~0.7), we can determine the 
contribution to the observed hadronic cross-section from the T(2S) and T(3S) resonances 
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FIG. 8. Mass recoiling against two charged particles, assumed to be pions, in events consistent 
with the kinematics for: e+e" jT{3S/2S)]T{3S/2S) T(lS)^+7r-; T(15) ^ /+/". Two peaks 
are evident; the leftmost peak corresponds to T(3S)^ 7r+7r^T(lS) transitions, the rightmost peak 
corresponds to T(2S)— > 7r+7r~T(lS) transitions. The calculated recoil mass differs from the true 
T(1S) mass due to our neglecting the (undetected) radiated photon in the recoil mass calculation. 

directly, by simply measuring the event yields in the peaks shown in Figure ^, and correcting 
by branching fractions and efficiency. 

For the contribution from the T(1S) resonance, we proceed by assuming that the photon 
spectrum varies as: ~ -g-, and that the production of a given Upsilon resonance is 

proportional to its dielectron width Fee. This gives a fairly simple prediction for the ratios 
expected for a given Upsilon resonance, since ~ (10.52 — Mr) GeV. We would expect 
that the production cross-section for T7 in e~^e~ annihilation therefore varies as: T{e~^e~ — > 
T7) oc This allows us to infer an expected production cross-section for 7T(1S') based 
on our measurements for 7T(2S') and 7T(35) production. We compare our extrapolated 
cross-section for e~^e~ — > 7T(1S) through this procedure with theory in order to estimate the 
magnitude of this correction. Combining our data with the theoretical predictions of Teubner 
et al. [||, we determine that the 7T(1S), 7T(2S), and 7T(3S) events comprise (1.8±0.6)% of 
the observed hadronic cross-section, where the error includes the uncertainties in the Upsilon 
decay branching fractions and detection efficiencies as well as the deviations between the 
theoretical and measured values. 

3. Two-photon collisions, which produce hadrons in the final state via e^e~ e~^e~'y'y 
e^e~-|-hadrons, are determined by running final-state specific 77 collision Monte Carlo events, 
and also by determining the magnitude of possible excesses in the EmsiUe vs. Ptransverse plane 
for data over qq Monte Carlo. These are determined to comprise (0.8ib0.4)% of our total 
hadronic event sample. 
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4. Beam-wall, beam-gas, and cosmic ray events are expected to have a flat event vertex distri- 
bution in the interval l^^tel <10 cm; such events are subtracted using a sample having a 
vertex in the interval 5.5cm < \zyrtx\ < 10cm, and extrapolating into the 'good' acceptance 
region <5.5 cm). These backgrounds are determined to comprise ~ (0.2 it 0.1)% of 
our hadronic sample. 

5. Remaining QED backgrounds producing 2 or more electrons or muons in the final state are 
assessed using a high-statistics sample of Monte Carlo events (to order in aQEo), and 
found to be < 0.1% of the sample passing the above hadronic event selection requirements. 

Summing these estimates results in a net background fraction / = (4.1 it 0.7)%. We note that, 
as this error is assessed partly by examining the difference between Monte Carlo hadronic event 
simulations and our data, this error also includes Monte Carlo modeling errors. 



IV. EFFICIENCIES AND RADIATIVE CORRECTIONS 

The computation of R is peformed with 

^ehad(l + '5)o-o^ ' 

where Nhad is the number of events classified as hadronic, / is the fraction of selected events 
attributable to all background processes, ehad is the efficiency for triggering and selection of events, 
6 is the fractional increase in hadronic cross section due to electromagnetic radiative corrections to 
that cross section, is the point cross section for muon pair production {86Mnb/E^^{GeV'^)), 
and C is the measured integrated luminosity. The luminosity is determined from wide angle e^e~, 
77, and iJ,'^fi~ final states and is known to ±1% For the data analyzed here, the integrated 
luminosity C is equal to (1.521 it 0.015) fb^^. 

To calculate R, we must therefore evaluate Eqn. ^. If the initial state radiation corrections 
were known precisely, we would be able to calculate the denominator term e(l -|- 6) with very 
good precision. However, since the uncertainties become very large as the center-of-mass energy 
approaches cc threshold {^/s ~4 GeV), the preferred procedure is to choose some explicit cut-off in 
initial state radiation (ISR) photon energy that makes us as insensitive as possible to the corrections 
in this high ISR photon energy/low hadronic recoil mass region. We therefore purposely design 
our selection criteria so that our efficiency for events with highly energetic ISR photons approaches 
zero. By choosing cuts that drive e to zero beyond some kinematic point, we ensure that the 
product e X (1 -|- (5) is insensitive to whatever value of 6 may be prescribed by theory beyond our 
cut. Thus, although there is a large uncertainty in the magnitude of the initial state radiation 
correction for large values of radiated photon momentum, we have minimized our sensitivity to 
this theoretical uncertainty. Figure displays our efficiency for an e'''e~ — > 'yqq event to pass 
our hadronic event criteria as a function of the scaled photon energy = „ ^ . We note that 
for >0.75 (corresponding to a qq recoil mass of Mrecoii < 5.25 GeV/c ), our integrated event- 
finding efficiency €had <1%- For x^ >0.75, we have therefore minimized our sensitivity to modeling 
uncertainties in this kinematic regime - increasing the initial state radiation contribution to this 
region results in a compensating loss of overall efficiency such that the product of e(l -|- 5) remains 
relatively constant. For our event selection criteria, we thus select the value of e(l -|- (5)(rE™"^ = 
0.75) = 0.90ib0.01, where the error reflects the systematic uncertainty in the radiative corrections. 
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FIG. 9. Efficiency for reconstructing e^e^ 'jqq event, as a function of the scaled photon 
momentum x-y = E^/Eheam- 

After subtracting ah backgrounds, dividing by the total luminosity, and normalizing to the 
mu-pair point cross-section, we obtain a value of i? = 3.56 it 0.01 (statistical error only). 

V. SYSTEMATIC ERRORS AND CONSISTENCY CHECKS 

We have checked our results in several ways. Backgrounds can be suppressed significantly by 
tightening the minimum charged track multiplicity to Nchrg >7, albeit at a loss of ~20% in the 
overall event-reconstruction efficiency. Imposition of such a cut leads to only a —0.4% change in the 
calculated value of R. Continuum data have been collected over 17 distinct periods from 1990-1996, 
covering many different trigger configurations and running conditions. We find a 0.3% rms variation 
between the various datasets used (the statistical error on R within each dataset is of order 0.1%). 
We can check contributions due to the narrow Upsilon resonances by calculating R using a small 
amount {5pb~^) of continuum data taken just below the T(2S) resonance, at £'6eam=4.995 GeV. 
We find that the value of R calculated using the T(2S) continuum agrees with that calculated using 
the T(4S) continuum to within one statistical error {Idstat)- Systematic errors are summarized in 
Table |. 

Systematic error summary 
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Source 


Error 


e X (1 + ,5) 


1% 


C 


1% 


Background Uncertainty/Hadronic Event Modeling Uncertainty 


0.7% 


Dataset-to-dataset variation 


0.3% 


Total 


1.8% 



TABLE I. Systematic errors in R analysis. 



VI. EXTRACTION OF as 

Using the expansion for R in powers of Os/vr given previously, with coefficients appropriate for 
this center of mass energy Q, we can evaluate the strong coupling constant. Using that expression, 
our value for R translates to 0^(10.52 GeV) = 0.20 ± 0.01 ± 0.06. 

The strong coupling constant can be written as a function of the basic QCD parameter A-p^, 
defined in the modified minimal subtraction scheme |^], 

47r / 261 Inix) 46? , , 1,, bobn 5,\ , , 



where 60 = (ll~2nj/3), /i is the energy scale, in GeV, at which as is being evaluated, 61 — ^ 



51-19ny^ 

5, = 2857 - ^ + X = /n(/i7A|^), and is the number of light quark flavors which 

participate in the process. To determine the value of q;s(90 GeV) implied by our measurement, 
we must evolve Ug across the discontinuity in ^jjg when the five-flavor threshold is crossed from 
the four-flavor regime. We do so using the next-to-next-to leading order (NNLO) prescription, as 
described in |^: a) We substitute 0^(10.52) into Eqn. (^) to determine a value for A-^ in the four- 
flavor continuum (obtaining A.j^{udcs)=4Q^ MeV). b) With that value of A-^g^, we can now again 
use Eqn. (^) to determine the value of as at the five-flavor threshold when the 6-quark pole mass 
(we use mb^poie='^-7 GeV) is crossed, and then use that value of as, as well as n/ = 5 in Eqn. (|3|) 
to determine ^j^s appropriate for the five-flavor continuum, c) Assuming that this value of A-^g^ 
is constant in the entire five-flavor energy region, we can now evolve Os up to the Z-pole, to obtain 
as{Mz) = 0.13 ± 0.005 ± 0.03, in good agreement with the world average as{Mz) = 0.118 ± 0.003 



VII. SUMMARY 

Near ^/s=10 GeV, R has been measured by many experiments, as shown in Table ||. The 
measurement of R described here is the most precise below the Z^. Theoretical uncertainties in QED 
radiative corrections (in the acceptance (ex (1 + 6)) and luminosity [^) contribute about the same 
amount to the systematic error as do backgrounds and efficiencies. Substantially improving this 
measurement will require progress on radiative corrections as well as on experimental techniques. 
Our R value is in good agreement with the previous world average, including a recent determination 



by the MD-1 Collaboration [18|. Our implied value of ag is in agreement with higher energy 



determinations of this quantity. 
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Experiment 



V^(GeV) 



R 



PLUTO II 



DASPII III 



DESY-Heidelberg |12| 



LENA m 



LENA HI 

cusB m 



CLEO 83 yj 



Crystal Ball |T 



ARGUS I? I 
MD-1 
Previous Expts., Weighted Average 
CLEO 97 (this work) 



9.4 
9.4 
9.4 
9.1-9.4 
7.4-9.4 
10.5 
10.5 
9.4 
9.36 
7.25-10.34 
w9.5 
10.5 



3.67 ±0.23 ±0.29 
3.37 ±0.16 ±0.28 
3.80 ±0.27 ±0.42 
3.34 ±0.09 ±0.18 
3.37 ±0.06 ±0.23 
3.54 ±0.05 ±0.40 
3.77 ±0.06 ±0.24 
3.48 ±0.04 ±0.16 
3.46 ±0.03 ±0.13 
3.58 ±0.02 ±0.14 

3.58 ±0.07 
3.56 ±0.01 ±0.07 



TABLE IL Summary of inclusive cross section measurements. 
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